Astronomy & Astrophysics manuscript no. paper 


February 2, 2008 


(DOI: will be inserted by hand later) 





Measuring Dark Matter Flows in Merging Clusters of Galaxies 

J. A. Rubino-Martm, C. Hernandez-Monteagudo, and T. A. Enfilin 
Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Str.l, Postfach 1317, 85741 Garching, Germany 



^1- 

o 
o 

(N 
X> 

PL, 

m 

(N 



> 

in 

(N 
in 

(N 

o 
o 

0,; 

6 ' 



Received date ; accepted date. 



Abstract. The Rees-Sciama effect produced in mergers of galaxy clusters is discussed, and an analytical ap- 
proximation to compute this effect from numerical simulations is given. Using this approximation and a novel 
toy model describing the physics of the merger, we characterize the spatial properties and symmetries of the 
Rees-Sciama signal. Based on these properties, we propose a method to extract the physical parameters of the 
merger, which relies on the computation of the quadrupole moment of the observed brightness distribution on the 
sky. The relationships between the quadrupole coefficients and the physical parameters of the merger (physical 
separation, projection angle on the sky and angular momentum) are discussed. Finally, we propose a method to 
co-add coherently the RS signals from a sample of cluster mergers, in order to achieve an statistical detection 
of the effect for those cases where individual signals are masked by the kinetic SZ effect, the primordial CMB 
components, and by observational noise. 
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1. Introduction 



The integrated Sachs- W olfe (lSW: ISachs and Wolfdll967t 
IHu and Sugivamalll994[) effect produces CMB tempera- 
ture fluctuations due to the accumulation of red- or blue- 
shift of photons travelling in time dependent gravitational 
poten t ials. Several r ecent s tudies llBoughn and Crittendenl 
120031 iNolta et all l2003t iFosalba and Gaztanagal |2003|) 
have claimed evidenc e for its detection in the WMAP data 
l|Bennett et al.ll2003|h The non-linear contribution to the 
ISW effect, in which the density contrast producing the 
gravitational potential is in its n on-linear regime, is usu- 
ally called t he Rees-Sciama (RS; iRees and Sciamalll968l 
ISelialJll99d) effect. 

Here we discuss the RS effect in the extremely non- 
linear regime of present day galaxy cluster mergers. The 
regime of moving single galaxy c l usters was already ex- 



amine d bv iBir kinshaw and Gull (1983), Aghanim^t^al 



1998 th iMolnar and Birkinshawl l|2000t l2003[h and lCoorav 



Galaxies: clusters: general - Galaxies: interactions - Gravitational lensing 



spectrally. However, both effects have distinct morpholo- 
gies, which in principle allow to discriminate them. 

The focus of this work is to examine methods to ex- 
tract the RS signature from merging pairs of galaxy clus- 
ters. The layout of the paper is as follows. In Sec. 2 we 
provide an approximation to compute the RS effect either 
from theoretical cluster models or numerical simulations. 
In addition, a brief description of the kSZ and the thermal 
Sunyaev-Zeldovich (tSZ) effects is also given, emphasiz- 
ing their phenomenological differences with the RS effect. 
In Sec. 3 we present and use an analytic toy model of a 
merger of two galaxy clusters to characterize the typical 
morphologies of these effects (RS, kSZ and tSZ). Based 
on its spatial properties, we present in Sec. 4 a method to 
extract both the RS signal and the physical parameters 
describing the system from a given merger. A method to 
stack the signal from a sample of mergers is also proposed. 
Conclusions are presented in Sec. 5. 



2002). However, galaxy clusters reach the largest veloci- 



ties, and therefore the strongest RS signal, during merg- 
ers in which two or more clusters are invoked. Due to the 
slow centre-of-mass velocity of the merging system, the 
RS signals of the merging subclusters can partly cancel 
or increase each other, depending on the merger geometry 
and the location of the line-of-sight (LOS). 

Furthermore, the spectral signature of the RS effect is 
a temperature change of the measured CMB photons, and 
therefore indistinguishabl e form the kinematic Su nyaev- 
Zeldovich (kSZ: ISunvaev and Zeldovichlll972lfl980|l effect 



2. Formalism 

2.1. Rees Sciama effect 

The CMB temperature change in a LO S direction n 
caused by the RS effect can be written as llCooravll2002T) 



ST RS (h) 



dr$(nr, r/c)T(n), 



(1) 



where <f>(r,t) is the time derivative of the gravitational 
potential and r/c the look-back time along the LOS. 
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In the following, we are dealing with well separated 
structures, galaxy clusters, which do not change signif- 
icantly during the time of a photon passage. Thus, we 
explicitly neglect the intrinsic variation of the cluster po- 
tential while the photon is passing through it. This allows 
an approximative treatment. 

In Newtonian approximation, the gravitational poten- 
tial depends on the total mass density g(r, t) via 



$(r,i) = -G fc 



r , Q(r',t) 
\r — r'\ 



(2) 



The time derivative of the gravitational potential therefore 
follows the change of the matter density, leading to 



$(r,t) = -G d 3 r' g(r',t)v(r',t) 



„/|3 



(3) 



by using the continuity Eq. j + V(«p) = 0, then the 
Gauss' theorem, and neglecting boundary terms in the in- 
tegration. Here, v is the mean velocity of matter, which is 
the density weighted mean of all matter species' velocities 
v = J2i v i Qil Qi where the i runs over all species like dark 
matter, gas, and galaxies, and g = J2i Qi- 

By inserting Eq. Elinto Eq.^we approximate all LOS 
to be parallel to the z-direction, we ignore the explicit 
dependency of g and v on time, and integrate out the 
LOS integral, yielding 



ST 



lis 



T 



M=-^ / d 3 r'g(r')v(r') 



(4) 



where vectors are split into a LOS-parallel and perpendic- 
ular part according to r = (r±,z), and where the Z-axis 
is taken pointing towards the observer. This approxima- 
tive description of the RS effect ignores small corrections 
due to non-parallel LOSs, gravitational lensing, the actual 
change of the matter distribution during the photon pas- 
sage, and the finite velocity of gravitational interaction. 

Nevertheless, its accuracy is sufficient to study the gen- 
eral difficulties in RS signal extraction from merging clus- 
ters of galaxies. Its simplicity allows easy inclusion into 
any code which calculates the kSZ effect from numerical 
simulations. 

One can introduce the LOS projected momentum 



P( r ±) = dz Q(r±,z)v(r±,z) 



in order to write 

STns 4G 
— (r±)=- — 



(5) 



(6) 



From this, one easily sees that the RS effect measures the 
convergence of the projected perpendicular momentum. 

2.2. Moving cluster of galaxies effect 

Eq. © can also be reproduced using a different ap- 
proach. The moving clust er of galaxies (MCG) effect 
(Birkins haw and Gulll Il98.^ was originally derived as a 



gravita tional lensing effect. How e ver, a nd as it was pointed 
out bv iMolnar and Birkinshaw! (|2000l) . this effect can be 
seen as an special type of RS effect, which is not caused by 
intrinsic variations of the gravitational field, but due to the 
movement of the cluster with respect to the rest frame of 
the CMB. Given that we have used exactly this assump- 
tion when deriving equation (@J, then our expression is 
completely expla ined as a gravitationa l len sing e ffect. 

Following iBirkin shaw arid Gulll ill flSajl and 
IMolnar and Birkinshaw! <|2000l) . the MCG effect for 
a cluster moving with velocity v, for small values of the 
lensing deflection angle, can be written as 



ST, 



MCG 



T 



(n) 



S(n) 



(7) 



where a is the angle between the cluster's peculiar ve- 
locity vector and its angular position vector, u v ± is an 
unit vector parallel to the projection of the velocity in the 
lens plane (v± = v smau v x), and 8(n) is the vector de- 
scribing the deflection angle . If we now use the Eq . for a 
geometrically-thin lens from ISchneider et al.l l|l992f) . 



S(r ± ) 



AG 



d 2 r'E(r' ± ) 



(8) 



where £ is the surface mass density, and we insert it in 
Eq. [7| we recover Eq. [B] in its limit of a constant velocity 
field. Since the MCG effect is linear, one can build a su- 
perposition of differently moving matter distributions, and 
therefore fully recover the more general Eq. [fj] It should 
be noted that Eq. [S] requires that the gravitational fields 
should be weak (which is true for clusters of galaxies), and 
that the matter distribution of the lens must be nearly 
stationary (i.e. w/c<C 1). 

2.3. Kinematic SZ effect 

Since the kSZ effect has the same spectral signature as 
the RS effect, it is crucial to include the latter into the 
considerations. Using similar planar and time-independent 
approximations the kSZ effect reads 



ST 



kSZ , 



T 



(r±) = dza T n e (rx,z)ve }Z (r±,z)/c, 



(9) 



where n e , and v e are the electron density and velocity, 
respectively. 

Since the electron and dark matter velocity fields are 
not too different, and the baryonic matter density follows 
roughly the dark matter density (g gas (r) s» /bg(r), with 
/b the cosmic baryon mass fraction), one can state that 



T 



fh - / \ 



(10) 



where m = m p {l + j A^ Hc )/ (1 + 1 -Xjie) is the gas mass per 
electron, Xn e = 0.24 the primordial He gas mass fraction, 
and m p the proton mass. 

Comparing Eq.[()]with Eci.llOlonc recognises properties 
which hopefully allow to separate the RS signal from the 
kSZ effect: 
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1. The RS effect is only sensitive to perpendicular mo- 
mentum, whereas the kSZ effect traces only parallel 
momentum. Thus, in a cluster merger which happens 
mainly within the sky plane the RS effect is enhanced 
relatively to the SZ effect. 

2. The kSZ effect is restricted to LOSs with baryonic mat- 
ter along it, whereas the RS effect reaches further out 
due to its gravitational origin. This can be used to 
enhance the relative signal strength. 

3. The RS has a more complex morphology compared to 
the kSZ since for a single moving mass, the RS has a 
dipole, whereas the kSZ a monopole morphology. For 
a pair of merging clusters with opposite momenta of 
the same strength this leads to a mostly quadrupole 
morphology of the RS effect and a dipole morphology 
of the kSZ effect. 

2.4. Thermal SZ effect 

Because of its different spectral characteristic the tSZ ef- 
fect can be separated from the the kSZ and RS effects. 
However, it has to be considered simultaneously not only 
because of possible residuals from incomplete removing of 
the tSZ effect from measurements, but also since the tSZ 
effect helps to localise the gas distribution and therefore 
to identify the regions which are expected to contain the 
strongest kSZ contamination. The usual tSZ y-parameter 
is given by 

y( r ±) = [dz<T T n c (r ± ,z) kTc ^ V± ' Z \ (If) 
J m c & 

where T e stands for the electron temperature, and mea- 
sures the LOS thermal gas energy content. 



Furthermore, hydrodynamical interactions of the gases of 
the two clusters are ignored, the gases are treated as non- 
interacting fluids. This coarse approximation is justified 
in so far that the kSZ and the RS effects both are only 
sensitive to momentum, and only in projection, but not 
to the thermodynamical state of the gas. The ignorance of 
shock waves exchanging the momentum of the gases there- 
fore does not strongly change the signal strength, only its 
exact spatial morphology. 

In contrast to that, the tSZ effect is very sensitive to 
the thermodynamical state of the gas and has therefore to 
be treated with a little bit more care. Again we do not try 
to model the exact spatial structure but try to get a good 
estimate of the signal strength. For that, we have to esti- 
mate the amount of shocked gas produced by the merger. 
We do this by measuring for each volume element the ex- 
cess kinetic energy of the gas of the two superposed clus- 
ters with respect to the volume's bulk kinetic energy and 
assuming this to be dissipated. The total thermal energy 
density of a volume element is therefore approximated to 
be 

etb(r) = 

(l d ° n ^ r ) fcr c,i(r-) + yn„(r) ( Vl (r)-v c (r))^j (12) 

where d c = 2 (1 + |A"h c )/(1 + ^Xh c ) is the num- 
ber of particles per electron, T e: i(r) (i £ {1, 2}) is 
the temperature of the gas in each cluster, and v c = 
SiGfi 2} v i n c,i/ Si n c,i ls the local mean gas velocity. 
From this, the electron pressure entering the tSZ effect 
estimate is easily obtained: 



3. Toy model 

3.1. Justification of approximations 

In order to be able to study the strength of the signal 
we construct a simplified toy model of a binary cluster 
merger. We grossly simplify the complicated gravitational 
and hydrodynamical process in order to have an analytic 
model. Nevertheless this model should provide us with 
maps of the RS, kSZ, and tSZ effects, which are suffi- 
ciently realistic for the purpose to study RS signatures 
by describing the typical strength and morphology of the 
effects. 

We build the merger model out of two cluster models, 
labelled with the numbers 1 and 2, and each des cribed 
by a NFW-DM matter profile l|Navarro et alJll997|) filled 
by isothermal hydrostatic gas. The two clusters approach 
each other. Each cluster velocity Vi (i € {1,2}) is cal- 
culated to be the velocity of a free-falling object in the 
gravitational potential of the other cluster. The total mo- 
mentum of the merging system is assumed to be negligible. 

Deformation of the clusters due to mutual tidal forces 
are ignored for simplicity, so that the DM and the gas 
of each cluster are streaming with the same velocity. 



(n e kT c )(r) = n etl (r)kT eil (r)+n ei2 (r)kT et2 {r) + 



m n cA (r)n ca (r) 
3d c rc e ,i(r) + n Cj2 (r) 



3.2. Simulations of the effect 



(v 1 -v 2 ) 2 (13) 



In this subsection, we present the equations outlined in the 
previous one for the simulation of a merger of two galaxy 
clusters, with a null total linear momentum. We assume a 
NFW dark matter profile for each cluster, 



Qdm(t) 



Qtt 



(r/r s )(l + r/r s y 



(14) 



The gas is considered isothermal, and distributed in hy- 
drostatic equilibrium, so we have IjMakino et alJ ll9981 



br s /r 



Qg{r) = Q 9 Q( 



1 



where b is a constant for each halo, given by 

_ 4:T:Gfj,m p g r 2 
ksT e 



(15) 



(16) 
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The value for g g o is determined using the following ex- 
pression for the gas mass fraction 

M gas (r v ) 



M to t(r v ) 



T- 



where the masses are evaluated at the virial radius 
(r„). We use here the va lues fi& = 0.044 and f2 OT = 
0.27 llBennett et all 12003^ . and we adopt T = 0.8 
l)Makino^r^dTTl998(i . It should be noted that when de- 
riving Eq. ^3 and when computing the total mass for Eq. 
IT71 we have neglected the gas contribution to the total 
density. 

For the computation of the RS effect, we use the 2-D 
projection of the NFW profile fEa. 1141) . which is given by 



When estimating from here the gravit ational energy asso- 
ciated to the two clusters, we follow ISarazinl (|2002) and 
use the approximation E grav w 1/2(</>iM 2 -|-</>2-Mi), where 
(17) is the gravitational potential of cluster i at the position 
of the centre of the other cluster. From this we calculate 
the relative velocity as 



(r±) = / dz Qdm{\J z 2 + r\) 




r± < r s 



r j_ > r s 



§ r s £>0j r± = r s 

A simulation of two merging clusters is then char- 
acterised by the set of parameters 8= {Mi, M 2 , d, i, J}, 
which correspond to the masses of the two clusters, the 
true physical separation, the projection angle of sky (i = 
90° is a face-on merger) , and the angular momentum of the 
system. Once these parameters are given, we can derive 
the remaining quantities describing the system. For defi- 
niteness, we use the following scaling relations. Both clus- 
ters are ta ken to be isothermal, with electron temperatures 
given by l(Eke et all Il99fit lAtrio-Barandela and Mucketl 
1999) 



T e = 7.75 



M 



2/3 



lO^h-iMr. 







(1 + z) keV 



Virial radii are derived as ijMolnar and Birkinshawll2000|) 

\ 1/3 



r v = 1.3/i _1 1 



M 







(1 + z)" 1 Mpc 



and the r« parame ter is computed as r c — 0.22r s 
l|Makino et alJ Il998j) . We have used a constant rat io 
r v /r c = 10 (self-similar core evolution, iKaiserl l|l98fil) 1. 
The electron density is derived as n e (r) = g g (r)(l + 
X)/(2m p ), where we fix the hydrogen abundance X to 
the primordial value (X = 0.76). 

Veloci ties are e stima ted following the formalism de- 
scribed in lSarazml l|2002() . although we modified it to use 
the gravitational energy associated to the NFW matter 
distribution described by Eq. El because it is more real- 
istic and better behaved at r = 0. For this latter profile, 
the gravitational potential is found to be 



4>(r) = —AttGqot 



3 ln(l+r/r.) 



2 

m 



E g rav{dQ^j -Egrav \d] 



1-'- 

d 2 



(19) 



(18) 



where d is the 3-dimensional separation of the clusters, 
do is the initia l distanc e from which clusters fall together 
(see Eq. 10 in ISa.ra.zinl l)2002lU . m = M 1 M 2 /(M 1 + M 2 ) 
is the reduced mass of the system, and b is the impact 
parameter of the merger. From here, the individual cluster 
velocities are calculated as v\ — M-ij (Mi+M^jv and v 2 = 
—Mi/(Mi+M2)v, so that the total momentum is zero and 
that v\ — V2 = v. 

3.3. Spatial properties/symmetries 

In this subsection, we use our toy model for two merging 
clusters to present two examples illustrating the spatial 
properties and symmetries of the effect. For definiteness, 
we choose the coordinate system with the Z-plane as the 
sky plane, with positive values pointing towards the ob- 
server. If we define ckcm as the angle between the direction 
joining the two centres of mass of the clusters (from Mi to 
M 2 ) and the X-axis, then we choose ctcM — (i.e. clusters 
on the X-axis, and Mi on the left of the image). This co- 
ordinate system will be used throughout the paper, unless 
otherwise stated. 

We first show in Fig.G]the result for a merger with 0= 
{1O 15 M , 5 x 10 14 M Q , 1 Mpc, 80°, 0}. This is a very opti- 
mistic case (in terms of the RS signal strength compared 
to the kSZ effect), where the merger is practically face-on. 

As pointed out in Sec. 2.3, the spatial pattern of the 
RS effect for a merging of two clusters shows a peculiar 
quadrupole structure. This quadrupole structure carries 
information about the velocity field of the clusters, and 
can exhibit, for example, if both clusters are in a pre- or 
post- merger state. 

This fact is used in the next section to extract the RS 
signal in a map where kSZ is also present, given that the 
kSZ has a dipole spatial structure. On the other hand, it 
can also be seen in the figure that the RS effect has a 
much more extended pattern than the kSZ, because of its 
gravitational nature. 

As a second example, in Fig. [21 we consider a more 
pessimistic case in which the velocity forms an an- 
gle of 60° with the line of sight. We incorporate also 
angular momentum. In this case, 8= {10 15 M Q , 5 x 
10 14 M Q , 1 Mpc, 60°, J}. Here a non-zero angular mo- 
mentum is implicitly defined by assigning a velocity direc- 
tion to the cluster Mi different from the direction joining 
the two cluster centres. The direction of Vi is specified 
by the spherical polar angles (8 V ,4> V ) = (60°, 60°), where 
9 V is the usual polar angle from the Z axis, and 4> v the 
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Fig. 1. RS (left image), kSZ (central image) and the sum of both effects (right image), for a pair of merging clusters 
with masses 10 15 Mq (left cluster) and 5 x 10 14 M© (right cluster), with zero angular momentum, when observed 
with projection angle 80°, and a physical separation of 1 Mpc. Contour lines in all plots correspond to (-64,-16,- 
4,-1,1,4,16,64) times the maximum of the map (dottcd-lincs for negative regions and solid lines for positive ones). 



usual azimuthal angle measured from the X-axis. In the 
case presented in Fig. [21 the RS effect is practically hidden 
by the kSZ, although it can be marginally seen at large 
distances. 

4. Signal extraction 

Once we have illustrated the morphology of the RS effect, 
we will study how to extract the RS information from 
observations of merging clusters. 

4.1. RS signature detection 

The peculiar structure of the RS effect in merging clusters 
can be quantified studying the dipole and quadrupole mo- 
ments of the observed intensity map. It is worth remarking 
that the quadrupole moment of the surface brightness is 
also used as a tool to extract the shear signal when study- 
ing weak lensing-induced d istortion of faint background 
galaxies ijKaiser et al.lll995j) . 

We propose the following recipe to proceed with future 
experiments aiming to measure this effect: 

1. Using a catalogue of cluster mergers, one identifies 
those clusters with a suitable projection angle (close 
to face-on mergers); 

2. Observations of the merger in frequencies far from the 
tSZ cross-over frequency (217 GHz, or 1.25 mm), or in 
X-rays, will show the spatial distribution of the gas. 
One of these images can be used to find the 'cen- 
tre of light' point, to which we refer the dipole and 
quadrupole moments 1 . In what follows, we will place 
the origin of our coordinate system at this point; 

3. Observations at 217 GHz (or performing a spectral 
component separation using observations at several 
frequencies) will provide us with a map where only 
kSZ and RS effects (as well as the CMB fluctua- 
tions) are present. In this map, we proceed to extract 

1 Note that, in general, the center of X-ray emission, the 
center of tSZ-light and the center of mass do not necessarily 
coincide in the same point. 



the RS component by computing the dipole (d) and 
quadrupole (Q) moments. 

We will adopt here the following definitions for the multi- 
pole moments of a given 2-dimensional map, M(x), with 
respect to the coordinate origin: 

d[M] = / ^-M(x)d 2 x (20) 
J \ x \ 

QiAM] = J 2XiXj ~l^ 2Sij M(x)<Fx (21) 

Weighting in this way, we assure that we give equal weights 
to equal areas. Thus, this is a robust approach in the sense 
that it should tolerate some level of noise (CMB or instru- 
mental) . 

For a 2D image, the Q-matrix only has 2 independent 
elements because is symmetric, and tr(Q)=0. Thus, the 
matrix is specified by quoting Q xx and Q xy . However, the 
quadrupole matrix can be rotated to its principal axis, so 
that a single number characterises the matrix (one of the 
eigenvalues), plus an angle specifying the orientation of 
the principal axis. It is easy to show that the two eigenval- 
ues for Q are ± ^Q xx + Q xy , which suggests to introduce 
the total, coordinate independent quadrupole moment as 

Q = \J Qxx "i" Qxy 

We will now illustrate how to extract information 
about the physical parameters describing the merger 
from the suggested recipe. For definiteness, we will as- 
sume here that we have two observational maps avail- 
able, one of the Compton y-parameter (y(x)), and an- 
other map of the brightness distribution on the sky at 
v = 217 GHz (B 2 n(x)). We will denote as Q tsz = 
Q[y(a;)] the quadrupole moment associated to the y- 
parameter map, and Q T = Q[_B2i7(a;)], the one associ- 
ated to the thermal temperature fluctuation map (Q T = 
Q RS + Q ksz + Q CMB ). We will start by considering the 
case of zero angular momentum, and neglecting the CMB 
contribution. These other cases will be discussed in the 
following two subsections. 
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Fig. 2. Same as Fig.^ but considering angular momentum, and with a larger projection angle. Here, i = 60°, and the 
velocity vector is denned by the spherical angles (60°, 60°). Contours have the same meaning as in FigQ] 



In Fig. |2| we show the dependence of the observed 
quadrupole with the projection angle. The coordinate sys- 
tem has been placed at the 'center-of-light' of the y- 
parameter map, and the maps have been rotated so the 
cluster centres in the tSZ map lie along the X-axis (so 
we only quote Q xx because Q xy = 0). Given that the 
weighting-factor in Eq. |2]is dimensionless, then Q T has 
units of flux, so in order to quote these values in mJy, 
we have assumed that the merging system is located at 
an angular distance of 300 Mpc. The map-size used for 
the computation of the quadrupole was 10.24 Mpc. This 
size has been used throughout the paper 2 . For illustra- 
tion purposes, the signal coming from the RS and the kSZ 
components has been computed separately, although the 
observable is the sum of the two quantities. Given that 
the RS signal has a more extended pattern than the kSZ, 
we have used a mask based on the tSZ map in order to en- 
hance the signal from the RS effect. In this example, this 
mask was obtained by tapering all pixels with a tempera- 
ture in the y-parameter map larger than 10% of the peak 
temperature. We can see that when using this taper, the 
quadrupole becomes totally dominated by the RS contri- 
bution. 

For illustration purposes, we also present in Fig. 0] 
the same plot but for the dipole modulus (d = \d\ = 

^Jd x + dy). We can see that this dipole signal is domi- 
nated by the kSZ contribution, and even when we mask 
90% of the tSZ emission (note that the dipole signal from 
kSZ is also reduced by approximately the same factor) the 
signals from kSZ and RS are comparable, but in this latter 
case only for low values of the projection angle (i < 50°). 
For (close to) face-on mergers, the kSZ contribution goes 
down because of the geometrical projection of the ve- 
locity (v z oc cosi) while the RS component is non zero 
(v± cx sini). Therefore, as expected, we conclude that the 
quadrupole moment is better suited for the extraction of 
the RS signal, while the dipole moment is more suitable 
for the kSZ (or tSZ) component (note that when consider- 

2 Note that one expects an small (logarithmic) dependence 
of the observed quadrupole with the map size. This can be 
easily seen for the case of two point-like clusters, where the 
asymptotic behaviour of the STrs signal goes as r~ 2 for large 
distances r from the center of mass of the system. 



ing only the inner part of the map, where 90% of the tSZ 
emission comes from, then the dipole is totally dominated 
by the kSZ component). 

In Fig. [3] we show the dependence of Q xx on the phys- 
ical separation, keeping fixed the other parameters. Here 
the separations are measured from Mi (the most massive 
cluster, i.e. the one with larger tSZ signal) to Ma, and the 
velocity for Mi is taken to be pointing towards the positive 
direction of X axis. Within this picture, a positive separa- 
tion (i.e. cluster Mi on the left) represents a pre-merger 
state, where the two clusters are approaching each other, 
while a negative separation (i.e. cluster Mi on the right) 
represents a post-merger state. Thus, we see that the sign 
of the quadrupole components carries information about 
the merger kinematics. We also stress here that tapering 
the central region of the map makes the quadrupole domi- 
nated by the RS signal. Finally, it is worth remarking that 
depending on the particular configuration of the merging 
system, the RS signal adds-up coherently or cancels in 
the computation of the quadrupole (or dipole), so we find 
that the outer region has a larger quadrupole-flux than 
the total map for some physical separations. 

4.2. Characterising the Angular Momentum 

The quadrupole pattern of the RS effect associated to 
merging clusters can also be used to learn about the angu- 
lar momentum of the system, at least in the plane of the 
sky: it is clear from Figs.^[21that the direction of cluster 
velocities conditions the nature of the quadrupole pattern 
imprinted by the RS effect. 

This is the motivation of the following approach, which 
is entirely based in the computation of the moments of the 
kSZ, tSZ and RS components, and hence is well defined 
for any morphology of the system: 

— First, one computes the quadrupole moment Q tsz (it 
would be valid to compute the direction associated to 
the dipole of the kSZ map, however, we shall ignore it 
due to its dependence on the inclination angle). From 
this moment, one computes the following angle, cor- 
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Fig. 3. Dependence of the measured quadrupole (Q xx ) on 
the projection angle i. These curves have been derived for 
the case of Mi = 1O 15 M (left cluster), M 2 = 5 x 1O 14 M 
(right cluster), separation of 1 Mpc, </>„ = 0° and the clus- 
ters placed along the X-axis (so Q = |Q ra |) and we use 
i = V (i.e. zero angular momentum). We performed the 
computation of the quadrupole in the whole image (quoted 
as 'all'), and masking the region that in the y-parameter 
map has a signal above the 10% of the peak intensity 
(quoted as 'outer'). We can see that when doing this, the 
quadrupole signal is dominated by the RS component. The 
flux is given in mJy, at the cross-over frequency of the tSZ 
effect {v = 217 GHz). The merging system is assumed to 
be at 300 Mpc. 

responding to one of the principal directions of the 
quadrupole: 



ctqtsz = arctan 



ftSZ 



ISZ 



+ {Q% s x z f 



.(22) 



Second, one computes the quadrupole moment Q T , 
which should be dominated by the RS effect due to the 
tapering of the central region. It is then convenient to 
take the angle related to the positive eigenvalue: 



Qqt = arctan 




(23) 



Taking as the origin of our coordinate frame the centre of 
tSZ-light of the system, the first angle should be related 
to the angle of the direction joining the centres of mass of 
each cluster (cccm)- Indeed, we find that aqtsz ~ acM- 
In a similar way, the angle chqt is related to the direction 
of the individual velocities of the clusters via Qqt ~ <Au/2. 
Thus, we can use the two observables given in Eg. 1221 and 
1231 to directly estimate the two physical angles olgm and 

This is illustrated in Fig. where we use the cluster 
parameters { Ml = 10 15 M Q , M 2 = 5 x 10 14 M©, i = 
9 V = 60°, d = 1 Mpc }, and we vary the azimuthal angle 



Fig. 4 

(d- 



Dependence of the measured dipole modulus 



idl + dy) with the projection angle i. We consider 

exactly the same case and notation as in Fig. [3| but here 
we also show the dipole signal from the region where the 
90% of the tSZ emission comes from (quoted as 'inner'), 
as well as the complementary 'outer' region. 





100 




50 













-50 




-100 




-2-10 1 2 

Separation (Mpc) 

Fig. 5. Dependence of the measured quadrupole (Q X x) 
on the separation between clusters. We consider the case 
Mi = 10 15 M Q , M 2 = 5 x 10 14 M Q , 4> v = 0° (so Q yy = 
0) and we use i = 6 V = 80° (zero angular momentum). 
Notation and units are the same as in Fig.|3l Separations 
are measured from the most massive to the less massive 
cluster, so in this picture a positive separation represents 
a pre-merger state, where the two clusters are approaching 
each other, while a negative separation represents a post- 
merger state. 

<fi v from to 2-7T radians. The quadrupole induced by the 
kSZ component perturbs our estimate of <p v from 2aqx 
as it is shown in the figure. Note that the value of <p v at 
which we have the maximum departure from the real value 
is not constant, but depends on the relative amplitude 
ratio of the kSZ and RS quadrupoles, i.e. depends on the 
other physical parameters of the merger. Finally, we again 
obtain better results for the case in which we mask out the 
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Fig. 6. Recovery of <f) v — acM from 2«qt — aqtsz . Thick 
solid line excludes that area containing 90% of the tSZ 
emission, whereas thin solid line considers the entire map. 
The dotted line gives the expected scaling 2aqx — aqtsz = 
<p v — ctcM- A projection angle of i = 60° was used in the 
calculation. 



clusters, it can be shown that Q^f = J z f{x\,X2), with 
x\,x% the coordinates of the cluster centres, and J z the 
projected component of the angular momentum along the 
Z-axis (see Appendix). Therefore, one can make use of 
the transformations of the components of the quadrupole 
under a rotation of the coordinate system by an angle a, 

Qxx = cos(2a) Q xx - sin(2a) Q xy 

Q xy = sin(2a) Q xx + cos(2a) Q xy . (24) 

in order to relate J z with the observable quantities 
(Q xx ,Q xy ) computed in an arbitrarily oriented reference 
frame. Using the observable aqtsz to estimate the orien- 
tation of the axis joining the cluster mass centres (cxcm), 
then the quadrupole components in the reference system 
aligned with the mass centre axis can be computed from 
Eq. I24I using a — aqtsz . As we have seen above, for high 
values of the projection angle i, we have that Q^ y w Q X y , 
and the latter is proportional to the magnitude of J z . In 
Fig. we illustrate this dependence of Q X y versus the Z- 
projection of the angular momentum, for the same physi- 
cal parameters as in Fig. H3 



0.10 



0.01 









all 




outer 


I 





\JJ ( 1 75 gr x cttl / s ) 



Fig. 7. Dependence of \Q xy \ with the absolute value of 
the projected component of the angular momentum, | J z \. 
We use the same physical parameters of the previous fig- 
ure, and computed for each value of <f> v the projected 
Z-component of the angular momentum (note that, for 
<j> v — 0° and <j> v = 180°, J = 0). Under the symmetry of 
our system (clusters along the X-axis), the kSZ residual 
affects mostly the Q xx values, leaving the crossed term 
Q xy practically unaltered. 

area where 90% of the tSZ emission comes from, even for 
cases with low values of the projection angle. 

With respect to the amplitude of the angular momen- 
tum, one can see that it is closely related to the cross 
component of the Q RS tensor, provided the X-axis of 
our coordinate system is located along the axis joining 
the cluster centres, and its centre is the system centre of 
mass. In this reference frame, a merger with J = shows 
no Q X y component. Indeed, for the case of two point-like 



4.3. Stacking of clusters 

The RS effect produces only a small temperature fluc- 
tuation making it very hard to measure it at individual 
clusters. It is therefore worth to develop ideas on how to 
co-add the signals from a large sample of clusters leading 
to a statistical RS effect detection. 

The simplest approach could be to co-add the 
quadrupole strength Q, which is always positive. However, 
there are two problems with this. First, the intrinsic CMB 
temperature fluctuations contribute also to Q, although 
their contribution can relatively accurately be estimated 
by measuring at non-cluster positions. Second, the kSZ 
effect imprints some quadrupole structure, which can not 
be completely masked away. Estimates of its strength are 
not possible from the observational data alone and there- 
fore would rely on numerical simulations of galaxy cluster 
merger in a cosmological setting. 

We therefore propose a different approach to stack 
the signal from several clusters. First, the cluster should 
be rotated in a way that its major axis (of the tSZ ef- 
fect or X-ray image) is aligned with the X-axis. Then, 
the quadrupole tensor of the temperature fluctuations 
qT = qCMB + qRS + Qksz ig ca i cu i a ted in this coordi- 
nate system and its components are co-added for all the 
clusters in the sample. We expect that for a sufficiently 
large number of clusters all contributions to Q T cancel 
out due to symmetry properties except the ones from the 
RS effect. 

The components of Q CMB have both signs with 
equal probability, since the CMB fluctuations are not 
correlated with the foreground cluster and the signs of 
quadrupole components are not rotational invariant (they 
reverse after a 90° rotation). For the same reason, we ex- 
pect that the lensing-induced cluster signal in the CMB 
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l)Seliak and Zaldarriagal2 000) will also cancel out, because 
the quadrupole moment is a linear function of the temper- 
ature, for which we expect both signs with equal proba- 
bility. 

Also the components of Q ksz have both signs with 
equal probability: mirroring a cluster through the X-Y 
plane going through its centre of mass reverses the sign of 
the kSZ map, and due to linearity also the sign of all Q ksz 
components. Since a cluster and its mirrored counterpart 
are equally likely in the cluster sample, their contributions 
to the co-added quadrupole cancel each other statistically. 

The components of Q RS have both signs with different 
probabilities. Since galaxy clusters are growing, one should 
find more clusters with converging velocity fields than with 
diverging ones. Diverging velocity fields should be present 
shortly after core passage of the dark matter clumps of 
two merging clusters. However, dynamical effects due to 
gravitational interactions (violent relaxation) should con- 
vert rapidly the organised bulk motion of the dark matter 
particles into an undirectional motion. This reduces the 
statistical contribution of diverging matter flows to the 
co-added Q T compare to the contribution of the converg- 
ing matter flows of early stage mergers. 

For this reasons, only a contribution from the RS effect 
due to the growing large-scale structure should survive 
statistically. However, there is one possible contamination 
to the described measurement scheme: any residual tSZ 
contamination to the temperature fluctuations is expected 
to be co-aligned with the major axis, and therefore to add 
coherently to Q T , unless they appear with both signs with 
equal probability. This point could be achieved by taking 
relativistic corrections to the tSZ effect into account. 

Therefore, further theoretical investigations are re- 
quired before this stacking method can be fully be estab- 
lished. This issue will be addressed in a follow-up paper. 

5. Conclusions 

In this paper we have presented a formalism that can be 
easily incorporated to N-body simulation codes in order 
to predict the Rees-Sciama effect in a merging systems of 
clusters of galaxies. For the typical range of velocities in 
those systems {v/c<€. 1), the obtained expression (Eq. 0} 
can also be seen as a gravitational lensing effect produced 
by a moving lens. 

Using simple modelling for the cluster merger, we have 
illustrated the morphology and symmetries of this effect, 
and we have developed a method to extract the signal, 
which should be applicable to realistic maps. This method 
is based on the computation of the (weighted) dipole and 
quadrupole moments of the brightness distribution. In 
particular, it has been shown how the quadrupole moment 
is related to the kinematic properties of the merger, so we 
can extract information about the dynamical state of the 
system: pre-merger or post-merger, and the magnitude of 
the angular momentum. 

Since we expect that in the near future the observa- 
tion of a single cluster merger will be extremely difficult, 



given the weak signal strength, we have proposed a simple 
method of stacking the signal from a large number of clus- 
ters in order to extract their RS signature statistically. The 
procedure is straight-forward: for a sample of clusters (e.g. 
a complete, or a merging cluster sample) a coordinate sys- 
tem is attached to the center of light (X-ray or tSZ map) 
so that the X-axis is aligned with the major elongation of 
the gas. Then the quadrupole moments of the CMB tem- 
perature fluctuations (outside the gas region, defined by 
the tSZ effect) are calculated for each cluster. Finally the 
individual quadrupole moment components are co-added 
for the sample. Intrinsic CMB fluctuations and kSZ ef- 
fect residuals should cancel out statistically, leaving only 
a signature of the RS effect. A detailed investigation of 
this stacking method will be carried out in a follow-up 
paper. 
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Appendix A: Relationship between the quadrupole 
moment and the angular momentum 

In this section, we prove that for two point-like clusters, 
Qxy is directly proportional to the component of the an- 
gular momentum projected along the line of sight, pro- 
vided the coordinate origin is located at the mass centre 
of the system and X-axis is aligned to the axis joining the 
clusters. Thus, the z=0 plane is the image plane. 

Let us assume that the coordinates of clusters are given 
by Xj = (xi,0), for i — 1,2. Their linear momenta will 
be labelled as = MjVj, i = 1,2 and we set the total 
momentum of the system to zero, so Pi = — P2. Also note 
that from the choice of the coordinate system, then Xi 
and X2 are related. The xy component of the RS-induced 
quadrupole is then given by: 



RS = 



-AG 



2 , 2xy 
ax - 



Px(x')' 



(x-x') 2 ' 



(A.l) 



with px(x') the projection on the z=0 plane of the LOS 
projected momentum density, given in Eq. 0] For our two 
point masses, we can take p_i_(x') = Pi<5d(x' — xi) + 
P2<5d(x' — X2), Sd denoting Dirac delta. It yields 



-8G I 2 xy 

a x — Pi 



x — Xl 
|x-xi| ; 



x - x 2 
|x-x 2 P 



.(A.2) 



The scalar product in this equation can be decomposed in 
two terms: 



■)RS 
Ixy 



-8G 



x — X\ 



(x — xi) 2 + y 2 



X ~ x 2 



(x — x 2 ) 2 + y 1 
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y 



(x - xi) 2 + y 2 

y 

(x - x 2 ) 2 + y A 



(A.3) 



The term proportional to P\. x vanishes when integrating 
over y, because the integrand is an odd function on y. The 
remaining term can easily be rewritten as: 



Q 



RS 

xy 



J 2 



(-8G) 



x y 
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(xi - x + x 2 - x) 



[(x - x x ) 2 + y 2 ][{x - x 2 ) 2 + y 2 



(A.4) 



with J z = —P\. v {x 2 — x\) the LOS projected component 
of the angular momentum. Therefore, 



Qxy = Jz f(xi,x 2 ), 



(A.5) 



with f(xi,x 2 ) a function depending exclusively on the ge- 
ometry of the system. 
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